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Abstract 

The recursion relations of hierarchical models are studied and contrasted with functional 
renormalisation group equations in corresponding approximations. The formalisms are 
compared quantitatively for the Ising universality class, where the spectrum of universal 
eigenvalues at criticality is studied. A significant correlation amongst scaling exponents is 
pointed out and analysed in view of an underlying optimisation. Continuous functional flows 
are provided which match with high accuracy all known scaling exponents from Dyson's hi- 
erarchical model for discrete block-spin transformations. Implications of the results are 
discussed. 
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1. Introduction 

Renormalisation group methods [1], and in particular Wilson's renormalisation group 
[2], play an important role in the study of physical systems at strong coupling and/or 
large correlations lengths. Differential implementations of Wilson's idea [3-10] rely on an 
appropriately-introduced momentum cutoff leading to flow equations for running couplings 
and A/^-point functions, which can be studied with a large variety of analytical and numerical 
methods. Numerical stability and reliability in the results is ensured through powerful 
control and optimisation techniques [11-14]. A different implementation of Wilson's idea is 
realised in hierarchical models of lattice scalar theories [15-19]. Hierarchical renormalisation 
group transformations are often discrete rather than continuous. Here, sophisticated 
numerical methods have been developed to extract the relevant physics, most notably for 
high-accuracy studies of scaling exponents for scalar models at criticality [20-22] and related 
theories (see [19] and references therein). 

Given the close similarity of the underlying principles, it is natural to ask whether 
Wilsonian (functional) flows can be linked explicitly, and on a fundamental level, to hier- 
archical models. If so, this link would provide a number of benefits. It will make powerful 
functional and numerical methods available to the study of hierarchical models. Vice versa, 
the numerical tools for hierarchical models could be employed for functional flows in specific 
approximations. Furthermore, an explicit link may lead to a path integral representation of 
hierarchical models, allowing for systematic improvements beyond a standard kinetic term. 
Finally, well-developed optimisation techniques for functional flows could be taken over for 
hierarchical models as well. 

In the limit of continuous hierarchical block-spin transformations, an explicit link between 
Dyson's hierarchical model [15, 17] and the Wilson-Polchinski flow [23] in the local potential 
approximation has been established long ago by Felder [24]. Following a conjecture of [13], 
this link has been extended [25, 26] to include optimised versions [11, 12] of Wetterich's 
flow for the effective average action [27]. These interrelations have recently been backed-up 
by extensive numerical studies of critical potentials and scaling exponents to high accuracy 
from either formalism [28]. 

In this paper, we evaluate the more general case and ask whether hierarchical models for 
discrete block-spin transformations are linked to functional flows with continuous renormali- 
sation group transformations. We first contrast the basic setups for functional fiows (Sec. 2), 
background field flows (Sec. 3), and hierarchical models (Sec. 4). At a Wilson-Fisher fixed 
point, underlying similarities and differences are worked out and compared for the leading 
scaling exponent (Sec. 5). An extensive numerical study of the eigenvalue spectrum of the 
Ising universality class from functional flows is performed (Sec. 6). A strong correlation of 
scaling exponents is established and analysed (Sec. 7). It is shown that specific functional 
flows match the leading and subleading scaling exponent from Dyson's hierarchical models 
for discrete transformation parameter to high accuracy (Sec. 8). We close with a discussion 
of the results and further implications (Sec. 9). 
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2. Functional Flows 



Wilsonian (functional) flows integrate-out quantum fluctuations within a path integral 
representation of quantum field theory. In their simplest form, they are generated through 
a cutofF term quadratic in the field added to the Schwinger functional, where the (classical) 
action is replaced by S — S+AS^ and AS^ ~ J dq(l){q)R}.{c^)(t){—q). The infrared momentum 
cutoff Rk{q^) ensures that the propagation of small momentum modes q^ <C fc^ is suppressed, 
while the large momentum modes (f' ^ k'^ remain unaffected. Under an infinitesimal change 
in the Wilsonian (infrared) cutoff scale fc, the effective action changes according to its 
functional flow, which reads (t = In k) 

dtVk = \tt {vf + Rk) dtRk (1) 

in the form put forward by Wetterich [27] . The trace denotes a momentum integration and a 
summation over fields. The factor dfRk in the integrand is peaked in the vicinity of q^ k, k'^. 
The cutoff function R^ obeys Rk{q^) — as k'^/q^ 0, Rk{q^) > as g^/A;^ — 0, and 
Rk{Q^) — >■ oo as A; — > A, and can be chosen freely elsewise, e.g. [11]. It ensures that the fiow 
is well-defined, thereby interpolating between an initial action at A; = A in the ultraviolet 
(UV) and the full quantum effective action T = Tk=o in the infrared A; ^ 0. 

In addition to providing a momentum cutoff, the function Rk{q^) also controls the 
stability and convergence of subsequent expansions [11-13, 29]. Therefore, it is possible to 
identify optimised momentum cutoffs - within given systematic expansions - which improve 
the physical result [11, 12, 30]. The construction of optimised cutoffs [11-14] is central to 
extract reliable results also in more complex theories including e.g. QCD [31], quantum 
gravity [32], thermal physics [5, 6, 33] and critical phenomena [28, 30, 34]. 

Below, we are interested in 3d scalar theories at criticality, where we can sent the ultraviolet 

scale A ^ oo. To leading order in the derivative expansion, the effective action reads F^. = 
/ d^x[\dfj^(f)dfj_(f) + Uk{p)] and p = ^((P'. Introducing r{y) = Rk{q^)/Q^ with y = q^/k"^, we find 

dtu = -3u + pu + dy (2) 
Jq y{l + r) + u' + 2pu" 

with u{p) = Uk{p)/k^ and p = p/k. An irrelevant constant originating from the angular 
integration has been rescaled into the potential and the fields. For the optimal cutoff i?opt = 
(A2 - g2)^(^2 _ ^2) ^i^i^ ^^p^ = (1/y - 1)^(1 - y), the flow reads [12] 

dtu = — 3it + pu' + 7^-— — - (3) 

after an additional rescaling. This fiow is integrated analytically in the limit of a large 
number of scalar fields [35]. We note that the universal content of the flow (3) is equivalent 
to the Wilson-Polchinski flow in the local potential approximation [13, 25, 26]. 
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3. Background Field Flows 



A different form of the flow (1) is obtained for momentum cutoffs wliich depend addi- 
tionally on a background field (p. Background fields arc most commonly used for the study 
of gauge theories [36], see [32, 37, 38] for applications. They have also been employed for a 
path integral derivation of (generalised) proper-time flows [39, 40]. 

In the presence of background fields, the functional Tj^lcp] turns into a functional of both 
fields, Tk[(l),4>]- In order to maintain the one-loop exactness of (1), the momentum cutoff 
can only depend on the background field, but not on the propagating field. Following [39], 



[0,(/)] = S''5'^Tk/5(j)''5q)"^. We 



we introduce x = F^'^'^'^ [(/),(()] and x = x[(j) = (j)], where r^"'™'' 
chose momentum cutoffs of the form R].{q^) xr[x\^ which depend now on the background 
field. Here, the regulator cuts off both large momentum modes ^ A;^ and large field 
amplitudes with r[.^'°^ ^ /c^. The full advantage of background fields becomes visible once 
they are identified with the physical mean, leading to the functional Tj^[(f),^ = (f)\ ^ ^ki'i^]- 
The resulting flow is closed provided Fj^. '[4>] = F^. ' [</>,</>]. For scalars, this relation becomes 
exact in the infrared limit studied below (for gauge fields, see [38]). Using the momentum 
cutoffs [39] 



rpT,m[ 



x\ 



exp 



we are lead to the background field flow 



dtTk = Tr 



k"^ + x/m 



2Fi[m,m;in + 1; 



4) 



x(l-|-rpT,m) \k'^ + x/m 



k' 



dfX . 



(4) 



If the term ~ dtx on the right hand side is dropped - meaning that additional flow terms 
originating from the implicit scale dependence in the momentum cutoff are neglected over 
the leading term - the flow (4) reduces to the proper-time flow of Liao [41]. A general 
proper-time flow is a linear combination of the first term in (4) for various m [40]; see [42, 43] 
for applications. 

Next, we specialise to the proper-time approximation to leading order in the derivative 
expansion. The flow equation for the effective potential takes the very simple form 



dtU = + pu' + ^ -r-nr , 

^ (m + + 2pM")™-3/2 ' 



(5) 



where m parametcrises the momentum cutoff, and an irrelevant constant factor has been 
rescalcd into the potential and the fields. For ra G [1, |], the fiow (5) is mapped onto the 
flow (2) [39]. At m = |, the flow (5) is equivalent to (3), modulo a trivial rescaling. As a 
flnal remark, we note that this proper-time flow is also obtained from linear combinations 
of higher scale-derivatives of Callan-Szymanzik flows, without relying on background fields 
[40]. In this representation, the approximation leading to (5) consists in the neglection of 
higher order flow terms 



(2) 
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4. Hierarchical Models 

Several hierarchical models for an effective potential v{(p) of a lattice scalar field have 
been introduced in the literature [15-17] (see also [19]). The hierarchical transformation 
laws relate the potential v{(p) at momentum scale k/i with an average in field space over 
v{(p) at momentum scale k, where £ > 1 is the renormalisation group step parameter. We 
restrict ourselves to the three-dimensional case; the generalisation to arbitrary dimensions is 
straightforward. 

In Dyson's model [15, 17], the renormalisation group step k k/£ for the potential is 
expressed as 

/ + 00 
dCfieiOe-'"'"'''^'"''^^^^ (6) 
-oo 

The details of the averaging procedure are encoded in the measure factor i^eiOi the ^- 
dependence of the potential on the right-hand side of (6) , and in the choice for the decimation 
parameter i. As is evident from (6), a decimation parameter £ = 2^/^ - employed for most 
numerical studies [19-22] - corresponds to a volume decimation of £^ = 2 at each iteration. 
For Dyson's model, the measure is chosen as fieiO = (^r o-(^))-i/2 exp(-^Vo"(^)) [19], where 
we require a{£) > for ^ / 1, and a{l) = with a'{l) 7^ 0. A standard choice is a{£) = 2(£—l) 
[24]. By definition, (6) describes a flow towards the infrared for decimation parameters £ > 1. 
For £ ^ 1, the hierarchical transformation (6) becomes continuous and the measure factor 
turns into a (5-function — > 5{^). Performing —idi{6), which is equivalent to kdk{Q), 
we arrive at a differential flow equation for the effective potential [24] 

dtv = -Sv + ^ifv' - v" + {v'f , (7) 

where an irrelevant factor is rescaled into the fields and the potential; t = Ink. The 
interaction terms in (7) originate from the scale-derivative of the measure —£diiJL^{^), which 
reads \a'{l)5"{5) in the limit £ ^ 1. This highhghts the relevance of the measure factor 
in hierarchical models. Our normalisation corresponds to the choice cr'(l) = 4 to match 
with [28]. The limit (7) is independent of a, but at £ 7^ 1, we expect that scaling solutions 
and exponents from (6) depend on it. Eq. (7) is the well-known Wilson-Polchinski flow 
[2, 23]. We therefore conclude that the potential in (7) is related to the potential in (3) by a 
Legendre transformation [26, 28]. 

A different version of a hierarchical model has been introduced by Wilson [16]. Here, the 
recursion relation is written as 

-00 

In Wilson's original model, the ^-dependence of the measure is = exp(— ^^), where 

the normalisation factor is ^-independent [19]. The measure factor is different from the 
one in Dyson's model, because the Gaussian width is independent. If instead we employ 
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the measure of Dyson's model, the limit i ^ 1 can be performed analytically.^ Up to a trivial 
rescaling, we find 

dtv = -3v + \^v' - v" . (9) 

In contrast to the Wilson-Polchinski flow (7), the non-linear term (-y')^ is absent. This comes 
about because the integrand of (8) as opposed to the integrand of (6) is manifestly 
symmetric under ^ Numerical evaluations of (8) have been reported in [16, 44]. For 

other representations of hierarchical models we refer to [19] and references therein. 

5. Matching Hierarchical Models 

In order to match hierarchical models by functional flows, we have to detail the scheme 
dependences of physical observables in either formalism. In the functional RG framework, 
the fully integrated flow is independent of the momentum cutoff Rk{q^) chosen for the 
integration. Scheme dependences, which enter as a consequence of truncations of [(/>], 
have been discussed extensively in the literature [11-13, 29, 45]. Their origin is easily 
understood. Since the momentum cutoff R in (1) couples to all operators in the theory, 
the missing back-coupling of operators neglected in a given truncation can result in a 
spurious dependence of physical observables evaluated either from Fq [(?!>], or from a fixed 
point solution F*[0]. The scheme dependence is reduced by identifying those momentum 
cutoffs, which, in a given truncation, lead to an improved convergence and stability of the flow. 

In Fig. 1, we discuss the scheme dependence quantitatively for the leading order scaling 
exponent u at a fixed point of the 3d Ising universality class [29]. Within exact flows (2), 
the full i?fc-depcndcnce has been studied in [30] by evaluating the fixed points of (2) for 
general cutoffs (Fig. 1, first column). The main result is that the range of achievable values 
is bounded from above and from below. The upper bound is attained for Callan-Symanzik 
type flows with ~ A;^. The lower bound with u = z/opt is attained with the optimal flow 
(3), and hence equivalent to the Wilson-Polchinski flow. The sharp cutoff result is indicated 
for comparison. 

The proper-time flow (5) rests on an intrinsically different truncation, because implicit 
dependences on the background field have been neglected as well as higher order flow terms 
proportional to the flow of F*^^); see Sec. 3. Therefore scheme dependences are quantitatively 
different. In the approximation (4), the m-dependence of scaling exponents from (5) has 
been studied in [39, 42] (Fig. 1, second column). The range of values is again bounded from 
above by a Callan-Symanzik flow. The lower bound is achieved for m — >^ oo. We note that 
the range of values exceeds those achievable within (standard) exact flows. The lower bound 



^ The variance of the Gaussian measure in (8) can be changed by an explicit rescaling of the fields as 
(p — > ^l\fo for finite cr, see [19]. Rescaling also the integration variable ^ and denoting the 

potential in terms of the rcscalod fields again as v{ip), we obtain (8) with a rescaled measure ^ie((,) ~ 
Ne t7~^/^ exp(— ^^/o"). It agrees with the measure of Dyson's model for Ne = n~^^^ and a = a{£). I thank 
Y. Meurice for e-mail correspondence on this point. 
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Figure 1: (colour online) Comparison of scaling exponent v from different functional flows (RG) in 
the local potential approximation, and hierarchical models (HM). The solid lines indicate the range of 
values obtained in the literature. The dashed lines, if present, indicate that the underlying parameter 
space has not been exhausted. The horizontal lines from top to bottom indicate the results for 
Callan-Symanzik flows, the sharp cutoff flow, the Wilson-Polchinski (optimal) flow, and the mean 
field result. A non-linear rescaling of the i^-axis is introduced for display purposes only (see main 
text). Colour coding: exact RG (red), exact background field RG in the proper-time approximation 
(violet), Wilson's hierarchical model (blue) and Dyson's hierarchical model (light blue). 



may be overcome once the additional flow terms, neglected here, are taken into account [39]. 
This is indicated by the dashed line. 

Next we consider scheme dependences of hierarchical models. Based on their construction, 
we expect that physical observables depend on the averaging procedure, on the measure 
factor /x^, and on the decimation parameter L It has proven difficult to systematically 
include wave function renormalisations and higher order operators in hierarchical models, 
and it is therefore not known whether the scheme dependence vanishes upon higher order 
corrections [18, 19]. Still, the scheme dependence should give a reasonable estimate for 
the underlying error in the model assumptions, in particular in comparison with functional 
methods. 

The ^-dependence of Wilson's model (8), originally constructed for £ 2, has been studied 
in [44] in the range I S [2^^'^, 2] (Fig. 1, third column). The full line covers the range of values 
obtained in the literature, while the dashed lines indicate that the underlying parameter space 
has not been exhausted. Because of (9) being linear in the potential as opposed to (7), we 
expect a strong ^-dependence, possibly a discontinuity, in the limit 1^1. The slope Idi v{tj 
along the data points with I > 2^/^ is negative, meaning that increases for smaller £. 
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We stress that Wilson's HM has an overlap both with exact flows and proper-time flows. 
Therefore, it is possible to map u{i) of Wilson's model for certain decimation parameters £ 
onto iy{R) from functional flows with appropriately chosen R. On the other hand, for some 
decimation parameter i, Wilson's model can only be mapped onto proper-time flows but not 
on exact flows, while for some decimation parameter it cannot be mapped onto either of them. 

The ^-dependence of Dyson's model (6) is displayed in Fig. 1, fourth column. The full line 
connects the known results at ^ = 1 [28, 30], i = 2^/3 [20, 22] and ^ = 2 [20]. The dashed line 
towards larger values for v indicates that the parameter space £ > 1 has not been exhausted. 
We note that the ^-dependence is very weak, with a tiny slope in the range of ^-values 
covered. The important observation is that the slope £d(,v{l) is positive in the vicinity of 
Ik. 1 — 2, implying that z^(^) > for ^ > 1. Consequently, it is possible to map the scaling 
exponent u{£) at discrete block-spin transformation i > 1 onto ^{R) from functional flows for 
speciflc momentum cutoff i?, both within the standard exact flows and within proper-time 
flows. This supports the conjecture that Dyson's model can be mapped onto functional flows. 

6. Spectrum of Eigenvalues 

Whether the observations of the preceeding section can be promoted to a full map 
between the formalisms crucially depends on further observables including the subleading 
scaling exponents. Here and in the following section, we study the spectrum of universal 

eigenvalues (scaling exponents) form functional flows (2) to high accuracy. A fixed point 
solution 7^ const, of (2) is characterised by the universal eigenvalues of cigcnpcrturbations 
in its vicinity. We denote the ordered set of eigenvalues as 0{R) = {ijJi{R), i = 0, • • • , 00}, 
with uji < LOj for i < j.^ In addition to the leading exponent i^{R) = — I/wq, we study the 
first three subleading scaling exponents lo{R) = ui{R), 0^2 and 0^3 (i?) within the exact 
fiow (2) for various cutoffs and coarse graining parameters. 

For the numerical analysis, we introduce several classes of momentum cutoffs deflned 
through wp = b/{{b+l)y-l), rcxp = l/(expc/-l); Vraod = l/(exp[c(y + (6- l)y^)/6] - 1), 
with c = In 2; and ropt,n = b(l/y — 1)"9{1 — y). These cutoffs include the sharp cutoff 
(6 ^ 00) and asymptotically smooth Callan-Symanzik type cutoffs R^ ~ k'^ as limiting 
cases. The larger the parameter b, for each class, the 'sharper' the corresponding momentum 
cutoff. The cutoff ropt,n probes a two-dimensional parameter space in the vicinity of Topt to 
which it reduces for 6 = 1 and n = 1. For integer n, ropt,n is a C*^"+^) function. In addition, 
we consider the cutoffs r,nix = cxp[-6(^ - l/^/y)] and rmix,opt = exp[-i(?/^ - y"'')], which 
obey rjnix(l/y) = l/rinix(2/)- Note that wc have covered a large variety of qualitatively 
different momentum cutoffs including exponential, algebraic, power-law, sharp cutoffs and 
cutoffs with compact support. Except for ropt,n, all cutoffs are C(°°)-functions. We employ 
the numerical techniques developed in [28, 30]. 



^ In our conventions, the sole negative eigenvalue at the Wilson-Fisher fixed point is wo- 
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Figure 2: (colour online) Six two-dimensional projections of the four leading scaling exponents in the 
Ising universality class from the functional flow (2) for various cutoffs and coarse graining (approx- 
imately 10"^ data sets). Here, v = — I/cjq. The Wilson-Polchinski result from the optimal flow (3) 
(large black dot) corresponds to a local extremum for all scaling exponents. Data sets based on rpower, 
J-mix.opt, ■Tmix, Texp, ^niexp, ''opt,! and r^od\ data from rpT is included in Fig. 4. 
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Our results for the universal eigenvalues at criticality are displayed in Fig. 2 for the six 
two-dimensional projections of the four-dimensional subspace {—000,10, 102, 1^3} of observables. 
The plot contains roughly 10^ data points, the different classes of cutoffs are colour-coded. 
We focus on the relevant 10%-vicinity of the Wilson-Polchinski result with scaling exponents 
Oopt = 0{Ropt) from (3), indicated by a large black dot, see [28] for the high-accuracy 
numerical values. The central result of Fig. 2 is that scaling exponents are very strongly 
correlated. Despite having probed the space of observables by many qualitatively different 
momentum cutoffs, we find that only a small subset of values can actually be achieved. The 
correlations increase the closer the eigenvalues O move towards Oopt- In the immediate 
vicinity of the Wilson-Polchinski result, we only find a very narrow "throat" connecting 
observables 0{R) with Oopt- For the sub-leading scaling exponents tOi, the throat remains 
very narrow even further away from Oopt- This is seen most clearly in the correlation of L02 
with tus in Fig. 2e), as well as in the correlations of u with both 002 and in Fig. 2d) and 
f). In turn, for the leading exponents v, the throat opens up more rapidly once its value is 
further away from fopt, see Fig. 2a), b) and c). 

7. Correlations of Eigenvalues 

The strong correlation of scaling exponents is a structural fingerprint of Wilsonian flows 
(2). Since the Wilson-Polchinski result is distinguished in the space of scaling exponents, it 
is natural to normalise the data of Fig. 2 with respect to it. We introduce the distance of 
any pair of scaling exponents {x,y) from the optimal result (xoptj^opt) as 

p{x, V) = ^(xopt-ar)2 + (yopt-y)2 = 10-^''(-'^) (10) 

We have chosen a standard metric in the space of observables (other choices can be applied 
as well). In this representation, full agreement with the (optimal) Wilson-Polchinski result 
is achieved ioi p ^ and Np ^ 00. We also introduce the angles 

ip{x,y) = arctan(a;/y) . (11) 

The critical indices u and ooi, i > I, are positive numbers. Therefore, they can cover the 
range x/y £ [0, 00] and ip G [0, f], and p > for any pair of observables {x,y). In the 
subspace {v,ijo), the extremal values (z^opt, "^opt) have the polar coordinates (popt, Vopt)) where 
the angle reads (y?opt = 0.78066 • • • which is close to 7r/4 = 0.785398 - - - , and popt = 0. The 
radial distance from the origin is p(0, 0) = 0.923002 - - - . 

In the representation (10), we can study the close vicinity of the Wilson-Polchinski result. 
In Fig. 3, we display our data points as functions of the angles ip, and their distance from 
Popt in a semi-logarithmic basis. It is noteworthy that only a very narrow range of angles ip 
is actually achieved by the data, despite the fact that large fractions of the underlying space 
of momentum cutoffs is covered. Also, and in contrast to Fig. 2a), many data points are 
degenerate in the representation {Np,(p). A priori, the Wilson-Polchinski value (popt, V'opt) 
could have been approached along many different paths. Instead, we find that only a narrow 
range of {Np, (^)-values is achieved for arbitrary momentum cutoff. 
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Figure 3: (colour online) Distance Np of the pair of scaling exponents {v, oj) from the optimal Wilson- 
Polchinski values (fopt, t^opt) in the representation (10), (11). Only a narrow range of angles ipivjUj) 
in the vicinity of « 7r/4 is achieved by the data. Many data points are nearly degenerate. Data 
from functional flows in the local potential approximation (2); same data sets and colour coding as in 
Fig. 2, plus further high resolution data points from ropt,i in the close vicinity of (fopt, Wopt). Results 
from the sharp cutoff limit, the Callan-Symanzik type flow (with i?^ ^ k"^) and Dyson's hierarchical 
model (with I — 2^/"^) are also indicated (black dots). The Wilson-Polchinski (optimal flow) result 
corresponds io ip — (/Sopt and Np oo. 



This pattern is further highlighted in Fig. 4, where we have magnified the non-trivial 
range of data sets from functional flows (2). In addition, we have added data points from the 
background field fiow (5) using the cutoff rpx^m foi^ ttl <^ ^. It is remarkable that these data 
sets display the same pattern as the data from (2). Our results from this and the preceeding 
section are summarised as follows: 

Extremum. — The Wilson-Polchinski (optimal fiow) result in Fig. 2 corresponds to an 
extremum in the space of physical observables with |ti;j| < Iwj^optl for all obervables in the 
vicinity of Copt- The extremum is local, because the exponents approach = 2i — 1, i > 0, 
for very soft (Callan-Symanzik-type) momentum cutoffs [30]. For the eigenvalue products 
n"^Q(a;j/wj^opt)) the Wilson-Polchinski extremum is a global one. 

Uniqueness. — Our result indicates that the correlations of eigenvalues at the Wilson- 
Pochinski result are strongest, in the sense that any flow of the form (2) with the exponent 
z^(i?) = Vopt automatically also agrees with the Wilson-Polchinski result in all other 
observables 0{R) = Oopt- In general, for J^(i?) > z^opti this is clearly not the case. 
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Figure 4: (colour online) Magnification of Fig. 3 in the vicinity of (/? ~ 7r/4 where v/u> k, 1. The data 
points for the distance Np{iy,uj) as a function of v/lu remain highly degenerate. Same data sets and 
colour coding as in Fig. 3, plus additional high resolution data points from background field flows 
(5) using rpT.m with m < |. Results from the sharp cutoff limit, the Callan-Symanzik type flow 
(with Rk ~ fc^) and Dyson's hierarchical model (with £ = 2^/'^) are also indicated (black dots). The 
Wilson-Polchinski (optimal flow) result corresponds to t'opt/wopt — 0.9905692 • • • and Np oo. 



Redundancy. — The eigenvalue correlations are so strong that the first two scahng 
exponents i^(i?) and uj{R), for a given R, contain enough information to fix the remaining 
observables on the percent level or below. These 'dynamical' constraints point at a major 
redundancy of (2) with respect to the underlying momentum cutoffs R. A relevant parameter 
has been identified previously. The gap minj^>oy(l + r) for normalised cutoffs [11], when 
maximised, leads towards the Wilson-Polchinski result [12, 13, 30]. 

Optimised observables. — Previous reasonings in favour of an optimisation only invoked 
properties of the underlying flow (1), e.g. its convergence, locality, stability and boundedness, 
allowing for improved physical predictions. This has been exemplified quantitatively for 
the observable v{R) which obeys 1 > i^{R) > fopt [30], where the lower bound fopt is 
closest to the physical result [25]. Fig. 2 now shows that this pattern extends to subleading 
eigenvalues. This equally extends to asymmetric corrections-to-scaling [46]. Therefore, one 
may turn the original reasoning around and argue that - because of the extremum property 
of the observables 0{R) - an extremisation of the functional fiow along the lines discussed 
in [11, 12, 14], or similar, should naturally lead towards the values Oopt- Stated differently, 
Figs. 2-4 show that observables derived from (1) admit an optimisation. 
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method 



cutoff parameter 



hierarchical model Dyson = 1) 



0.649 561 773 880" 0.655 745 939 193« 



functional RG 



r-opt,n (n = 1, 6 = 1) 

^compact {p ^ 0) 

nnt {h ^ 1) 

r-pT,m (m = 5/2) 



0.649561773880" 
0.649561773880 
0.649561773880 
0.649561773880" 



0.655 745939193" 
0.655 745939193 
0.655 745939193 
0.655 745939193" 



Table 1: Matching scaling exponents v and u from continuous hierarchical transformations with 
functional flows. Results agree at least to the order 10~^^. Data from this work, and from a) [28]. 



Finally, we note that the data point from Dyson's hierarchical model with I = 2^^^ - 
as plotted in Figs. 3 and 4 - nicely fits into the set of data points covered by functional 
flows, extending the link observed in Sec. 5 beyond the leading exponent. This observation 
is addressed quantitatively in the following section. 



8. Matching Beyond the Leading Exponent 

To further substantiate our conjecture that hierarchical models could be mapped onto 
functional flows, we have to show quantitatively that results from hierarchical model are 
reproduced by specific functional flows. Here, we study the close vicinity of the Wilson- 
Polchinski result Oopt; where the correlations are strongest, see Figs. 2 and 3. We have to 
restrict our search to Dyson's hierarchical model, where high-accuracy data for the first 
sub-leading scaling exponent u is available. No subleading exponents have been computed 
for Wilson's model. 

For the numerical analysis, we introduce additional classes of momentum cutoffs R which 
contain the optimal flow (3) in some limit. In addition to the two-parameter family of cutoffs 
'"opt,™) we also study the cutoff rcompact = exp[— e~"'^/^/(6 — y)\6{h — y) for 6 > which 
is C^°°\^ and the cutoff r;nt = cxp(— — y)0{y — b) with b G [0, 1], which is effective for 
a finite interval of momenta G [bk'^,k'^]. In the limit 5^0 (6 — > 1), the corresponding 
flows are equivalent to (3). Hence, ropt,m ^compact and Tint parametrise substantially different 
classes of cutoffs. More generally, there are infinitely many cutoffs Rj- leading to scaling 
exponents identical with Oopt, and the examples provided above serve to illustrate this. 

At ^ = 1, Dyson's hierarchical transformation is continuous, and the scaling exponents 
are equivalent to those from the optimal flow (3) and the Wilson- Polchinski flow (7). In 
Tab. 1, we compare exponents from different functional flows. We confirm numerically, and 
with high accuracy, that the cutoffs ropt,m ''compact and Tint lead to the Wilson-Polchinski 
result for specific parameter values. 
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method 


cutoff 




parameter 


v 


Ul 


hierarchical model 


Dyson 


{i = 


2V3) 


0.649 570'' 


0.655 736* 






in = 


1,6= 1.048) 


0.649 570(9) 


0.655 736(6) 






(n = 


1,6 = 0.9545) 


0.649 570(9) 


0.655 736(9) 






in = 


1.135, 6=1) 


0.649 570(6) 


0.655 736(8) 




' opt,n 


(n - 


11 h — -t r)28~i 


649 570f6'l 


655 7'^6('8'l 


functional RG 


^compact 


{b = 


0.04775) 


0.649 570(9) 


0.655 736(9) 




Tint 


(6 = 


0.944) 


0.649 570(9) 


0.655 736(8) 




Tint 


(6 = 


0.9444) 


0.649 570(7) 


0.655 736(9) 






(to = 


= 2.499785) 


0.649 564(9) 


0.655 736(1) 




rpT,m 


(to = 


= 2.49944) 


0.649 570(1) 


0.655 720(6) 



Table 2: Matching scaling exponents v and lu from discrete hierarchical transformations with functional 
flows. Results agree to the order 10~^ for all cutoffs except the proper-time flow, which matches up 
to the order 10~^. Data from this work, and from b) [22]. 

At ^ = 2"^/^, Dyson's hierarchical transformation is discrete. The reference data reads 
'^DHM ~ 0.649570 and w^jj^ = 0.655736 [22].^ These values differ only at the order 10~^ from 
the optimal (Wilson-Polchinski) result, and are therefore sufficiently close to Oopt to confirm 
or refute the correlations observed in the previous section. Figs. 3 and 4 indicate that the 
result from Dyson's hierarchical model is fully matched by functional flows. Our numerical 
results are given in Tab. 2; brackets indicate that a digit is possibly affected by numerical 
errors. We have found several sets of parameter values, such that the scaling exponents 
agree with Dyson's model to order 10^^. More importantly, the momentum cutoffs arc quite 
different. Hence, our analysis also confirms the strong correlation of scaling exponents in 
the immediate vicinity of the Wilson-Polchinski result. Based on the eigenvalue correlations 
within functional flows, we conjecture that the subleading eigenvalues Ui with z > 2 of 
Dyson's model at ^ = 2^/^ also agree to the corresponding accuracy with the values implied 
through the functional flows in Tab. 2. 

Our results based on the proper-time flow (5) with rpT has also been given in Tab. 2. Full 

agreement is achieved either with the exponent z^dhm °^ the subleading exponent w^hm' but 
not with both of them. Once one of them is matched, the deviation in the other observable 
is of the order 10^^. The relevant parameters are m < |, the regime where (5) is mapped 
onto (2) [39]. Therefore, the values in Tab. 2 reflects well the range covered by standard 
Wilsonian flows (2). We expect that full agreement is achieved for proper-time flows which 
are linear combinations of (5) for different m, but we did not attempt to do so here. 



■* In [22], high-accuracy results at £ — l}^^ have been given for = 1v and ^ = vu) (and = 0) with 13 
significant digits. They imply (/dhm = 0.649570365 • • • and wdhm = 0.655736286 • • • . For the present study, 
only the first six figures are required. 
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In summary, we have provided numerical maps from several functional flows onto 
Dyson's model at a non-trivial i ^ 1, with an accuracy of the order lO^^. The set of 
achievable values for scaling exponents from functional flows in the close vicinity of the 
optimal result is just wide enough to accommodate for the data from Dyson's model. 
This is a non-trivial result, also showing that the ^-dependence of Dyson's model and 
the i?fe-dependence of functional flows are very intimately related. Based on our results 
for v and a; at ^ 7^ 1, and on continuity in i, we expect that this map extends to other 
universal quantities in the same approximation, analogous to the full map which is known 
for £ = 1. Data for further symmetric and asymmetric corrections-to-scaling exponents, 
once available, will allow for additional checks of this picture. Full equivalence is guaranteed 
as soon as an explicit link in the form i = £{Rk) or Rk{q^) = Rk{q'^,£) is furnished. For 
the local potential approximation, our results indicate that this map, if it exists, is not unique. 

9. Discussion and Conclusions 

Establishing equivalences between implementations of Wilson's renormalisation group as 
different as discrete hierarchical models of lattice scalar fields on one side and continuous 
functional flows on the other, allows for new views and insights on the respective formalisms 
and on the underlying physics. Previously, equivalences were known only in the limit where 
the hierarchical transformation becomes continuous. In this paper, based on similarities 
in the dependences related to the underlying coarse-graining, we have extended this link 
towards discrete hierarchical transformations. This correspondence shows that continuous 
RG flows (1) are sensitive to implicit discretisation effects via the momentum cutoff. 

Specifically, for the 3d Ising universality class, we have compared the formalisms on the 
level of scaling exponents. Their dependence on the step-size parameter i within Dyson's 
hierarchical model (6) is qualitatively and quantitatively similar to their dependence on 
the momentum cutoff R within functional flows (2). In cither case, scaling exponents 
are bounded by the Wilson-Polchinski values obtained for £ ^ 1 and R — > -Ropt- Once 
the hierarchical transformations are discrete, £ ^ I, slight variations in all known scaling 
exponents from Dyson's model are matched by functional flows with non-optimal momentum 
cutoffs R ^ Ropt- This is quite remarkable, particularily in view of the strong eigenvalue 
correlations found amongst functional flows. In this light, the optimisation of functional 
flows with R Ropt can now, alternatively, be viewed as the removal of discretisation effects, 
at least to leading order in a derivative expansion as studied here. It will be interesting to 
contrast these findings with the construction of improved or perfect actions on the lattice. 
More generally, it is conceivable that Dyson's model for arbitrary £ is mapped by functional 
flows on a fundamental level beyond the numerical map provided for i = 2^/^. An explicit 
map would be very welcome, also in view of linking hierarchical models to a path integral 
representation of the theory. In Wilson's hierarchical model (8), the range covered by the 
leading scaling exponent indicates that a partial map onto functional flows exists, though 
only for a restricted domain of ^-values. Interestingly, the overlap with background field 
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flows is even larger. Whether these maps extend beyond the leading exponent cannot be 
settled presently due to a lack of data for subleading exponents from Wilson's model. 

In addition, we found a distinct correlation of scaling exponents from functional flows 
(2). The eigenvalue spectrum, a fingerprint of the physics in a local potential approxima- 
tion, is severely constrained and achieves the Wilson-Polchinski values as an extremum. 
Furthermore, the full space of physical obscrvablcs is described by very few parameters only, 
instead of the infinitely many moments of the momentum cutoff. Resolving this redundancy 
should prove useful for studies of e.g. non-trivial momentum structures and higher orders in 
the derivative expansion. Finally this pattern amongst physical observables highlights the 
extremum property of cutoffs leading to the optimal Wilson-Polchinski result [11-14]. We 
expect that the intimate link between optimised flows on one side, and extremum points in 
the space of observables on the other, persists in more complex theories. This observation 
will prove useful for studies in QCD and quantum gravity, where an appropriate optimisation 
is even more important to extract the relevant physics. 
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